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Abstract 

The electromagnetic form factors are computed using eigenstates of linear mo- 
mentum for the nucleon. The latter is described in the framework of the chiral 
color-dielectric model, projecting the hedgehog ansatz on eigenstates of angular 
momentum and isospin. Form factors are well reproduced, with the exception of the 
magnetic one for the proton, up to q 2 ~ 0.5 GeV 2 . The effect of the removal of the 
spurious center-of-mass contributions shows up mainly in the electric form factor of 
the proton. A noticeable improvement is obtained with respect to the calculation 
without linear momentum projection. 



1 Introduction 



We report on a theoretical calculation of the electromagnetic form-factors of 
the nucleon in the space-like region, performed in the framework of an effective 
model — the chiral color-dielectric model (CDM) [|1],|2| — in which the nucleon 
is described as a chiral soliton. The model contains quark and meson degrees 
of freedom and a phenomenological scalar field which is responsible for quark 
confinement. 



In the previous calculations of the form factors in the framework of the chiral 
soliton models of the nucleon, as the linear sigma model , the Nambu- Jona- 
Lasinio model f§], the Skyrme model f| and the CDM ||, it was always 
assumed that the nucleon is at rest before and after the interaction with 
the virtual photon, the so-called static approximation. In the present work 
we overcome, at least in part, the technical difficulties associated with the 
computation of the form factors when the nucleon initial and final states are 
eigenfunctions of linear momentum, at least non relativistically. Our formalism 
is a generalization of the one presented in a work by Neuber et al. [0|§, where 
static properties of the nucleon have been computed in the framework of the 
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CDM. In their case therefore it was enough to build eigenstates of the angular 
momentum having zero linear momentum. 

The main technical problem in our calculation is due to the non-commutativity 
of the projectors on linear and on angular momentum. It will be shown that 
the problem can be solved by taking a suitable average on the direction of the 
transfered momentum, as it is suggested by the equation / dq [P q , Pjm] = 0. 
Morover it will be shown that the Fourier transform of the matrix elements of 
the electromagnetic current does not depend on the direction of the momentum 
transfered, if nucleon states are considered. Therefore the integration on q 
does not imply any approximation. The use of the effective commutativity of 
the two projectors simplifies the computation of the doubly projected form 
factors. These mathematical aspects are applicable in general to any quark- 
meson chiral soliton model. 

The numerical results shown in this paper refer to the electric and the magnetic 
proton and neutron form factors computed in a particular version of the chiral 
CDM. This is the so-called "single minimum" version which gives good results 
in quark matter calculations ||. The model contains two parameters that we 
adjust in order to reproduce the average A—N mass and the isoscalar nucleon 
radius. All other results are obtained without any further parameters' fitting. 

This paper is organized as follows. In Section 2 the electromagnetic form- 
factors are defined. In Section 3 we review the chiral color-dielectric model 
and the way model states representing a nucleon with definite momentum are 
obtained. Section 4 is devoted to the formalism to compute the electric and the 
magnetic form factors of the nucleon in the projected hedgehog state. Finally, 
in Section 5, the results are presented and discussed. 



2 The electromagnetic form factors of the nucleon 

Let l-ZVad^)) represent a nucleon state of mass Mm, with spin and isospin 
described by a. The standard definition of the electromagnetic form factors is 
given by: 

(N f (p' ll )\JLmN l (p ll )) = u f (p[ l ) 



They only depend on the modulus of the momentum transfer of the virtual 
photon q^, = p'^ — p^, being real functions of q^ = (q ) 2 — q 2 . We work in 
the Breit frame where the photon 4-momentum is g M = (0,q), i.e. the energy 
transfer is zero. Our results will be presented as a function of q = \q\. 
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In the Breit frame, Eq. ([!]) reads explicitly 



Wf) I -C(o)|iVi(-f )> = %(§) 



«i(-f)-(2) 



Normalizing the Dirac spinors so that uu = 1, 

( , 



IE + M. 



u{p) 



N 



2M 



N 



(3) 



the matrix element (§) can be worked out yielding 



(N f ^)\JL(0)\N t (~p 



- F2(q 2 ) 



(X/ IXi/^o + « 2^ <W 

{XflXt UM 2 N »° h 2M N 



(4) 



From Fi(q 2 ) and F 2 (q 2 ) it is usual to define the so-called Sachs form factors, 
Ge{q 2 ) and Gm{<1 2 ), where E and M stand for "electric" and "magnetic" 
respectively, which are expressed by 



G E (q 2 )=F 1 (q 2 



G M (q 2 )=F 1 (q 2 ) + F 2 (q 2 



(5) 
(6) 



Using these definitions and Eq. (|4]) one obtains explicit formulas for the electric 
and the magnetic form factors: 



G E (q 2 )(Xf \Xi) = W§) I J l(0)| Ni(-p 
.G M (q 2 



9 Uy -(Xf W x q\ Xi ) = (N f q) | J em (0)| JVi(-|)> . 



(7) 
(8) 



3 The projected chiral color-dielectric model 



In this work we shall use a chiral version of the CDM, whose Lagrangian reads 
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C = iqYdf,q + -q (a Q + i^ 5 r ■ tt) q + -{d^xf - V (x) 
X 2 

+ ^(d,a ) 2 + ^(d,n) 2 ~U(a ,7r). (9) 

This is a model with interacting quarks, chiral mesons a Q and 7r and also a 
chiral singlet scalar field % responsible for confinement. The potential U (a Q , n) 
in (|9|) for the chiral mesons is the so-called 'mexican-hat' potential 

A 2 2 
U (a Q , tt) = — (a 2 + tt 2 - u 2 ) + ca a + d. (10) 



The chiral symmetry SU(2) x SU(2) of L is explicitly broken by the small 
term ca Q in (|T0"D; the last term in the same expression is a constant fixed in 
order to have min U = 0. The parameters A and v are related to the sigma 
and the pion masses and to the pion decay constant: 



ml-ml_ 2 2 ml 



y= \ u ~ i u — U~ p • (11) 



For the x field potential we consider the quadratic form 

V(X) = \m 2 X 2 , (12) 



where M is the x mass. It is well known that the chiral CDM allows for soliton 
solutions in which the quarks are absolutely confined |l|,@J3l • m sucn solutions 
the x mean field is a decreasing function of the distance, approaching zero 
in the limit r — > oo. This generates a raising dynamical mass for the quarks 
and confines them. In previous works an exhaustive study of the model with 
a quartic (or 'double minimum') potential was carried out In this work 
we consider just the quadratic ('single minimum') potential for the confining 
field. We recently showed that for 'double minimum' potentials and for all 
sets of parameters fitting nucleon properties, the equation of state for quark 
matter turns out to be unrealistic. Indeed even at very low density the energy 
per baryon number for quark matter turns out to be smaller than that for 
nuclear matter. Using instead a quadratic potential for the x field a realistic 
equation of state is obtained. 

Altogether, the parameters of the model defined by @ are the pion and sigma 
masses (fixed at m n = 0.139 GeV and m a =l. 2 GeV), the pion decay constant 
(/„. = 0.093 GeV), and g and M, the quark-meson-x coupling constant and 
the x-mass, respectively. 

In order to obtain model states representing the nucleon we used the pro- 
cedure explained in great detail in Ref. J7| which, for the sake of complete- 
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ness, is sketched here. We consider three valence quarks in the hedgehog state 
\h) = (\u I) — \d ])), all occupying the same lowest positive energy s- 
orbital and surrounded by clouds of \i sigmas and pions, described by coher- 
ent states ( | n), for pions; | S), for sigmas; and | \)i f° r the confining field). 
The meson mean fields are the expectation values of the field operators in the 
corresponding coherent states. 

We can write the quark single particle states and the meson mean fields as 



= ~tr= [ ..iQ A )\h) (13) 



/4tt \iv{r)cr ■ r , 

Vo{r) — (S \cr\ £} — /tt = cr(r) — (14) 

7r(r) = (niTrin) = -Mr) (15) 
r 

X(r) = (x \x\ X) = X(r) , (16) 

where a is just the fluctuating part of the a Q field. Altogether, the hedgehog 
ansatz reads 

|^) = (cl) 3 |n)|s)|x), (17) 



where Cl creates a particle in the single quark state fljjp. 

Of course, solitons described by the hedgehog | iphh) cannot represent physical 
baryons because they are not eigenstates of angular momentum or isospin. 
In addition, ( |TT[ ) represents a localized object and therefore the translational 
symmetry of the model hamiltonian is also broken in such states. In particular 
they contain spurious centre-of-mass components which contribute to the en- 
ergy and to the other observables. However, a nucleon at rest can be obtained 
by applying the projector onto linear momentum q = together with the 
projector onto angular momentum-isospin. The linear momentum projector is 
given by JlO| 



where U(a) is the translation operator. It is well known that due to the sym- 
metry of the hedgehog it is enough to perform a single projection (e.g. onto 
spin), since this automatically projects onto the same value of isospin. The op- 
erator which projects out from the hedgehog a state with angular momentum 
J and isospin T = J is 

P J mm t = (- 1 ) J+MT ^^ J d 3 nvil_ MT (Q)R(Q), (19) 
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where Q = (a,/3, 7) are the three Euler angles, V J M K {VL) are the Wigner 
functions and R(Q) is the rotation operator. In the following we will consider 
M T = —M and use the shorthand notation P JM = P^-m- 

The radial functions in (|!3| )-(ll6|) are determined using an approximate variation- 
after-projection method firstly suggested by Leech and Birse |[L1|| . They are 
computed by minimizing the expectation value of the (normal-ordered) model 
hamiltonian in the model baryon state with quantum numbers J — T — | 
and linear momentum zero: 

P q =oPjM I *hh ) = I J, T, M, q = ) . (20) 



In the model there are two parameters, g and M, yet to be fixed. However, 
because of the smoothness of the x-held in a typical soliton solution and of 
the relative weakness of the chiral meson clouds, the relevant parameter turns 
out to be G = \fg~M. In the quark matter sector this is indeed the only free 
parameter of the model [12[]. If G is fixed to reproduce the isoscalar radius 



of the nucleon one obtains G = 0.2 GeV. For this parameter the nucleon- 
delta average mass is around 1.13 GeV (experimental value 1.085 GeV). It is 
interesting to observe that if g and M are changed, keeping G fixed, the static 
properties of the nucleon are essentially unchanged. For example, for G = 0.2 
GeV, changing the mass of the x field in the range 0.8-2.0 GeV, affects the 
results by less than 1%. 

In the present version of the CDM, the nucleon-delta mass splitting results 
only from the quark-pion interaction. Due to the weakness of the pionic field, 
the nucleon delta mass splitting obtained is too small. The experimental value 
of the splitting could be recovered if, in addition, a color-magnetic interaction 



(like in the MIT bag model or in the cloudy bag model) was considered |13 
We will come back to this point in the conclusions. 



4 Electromagnetic form factors in the projected hedgehog 



In order to compute the electromagnetic nucleon form factors one has to eval- 
uate matrix elements of the electromagnetic current operator. In the CDM 
the latter is given by 

JU*) =■ E Q + I) q a (x) + [tt(x) A d^(x)} 3 : . (21) 

As mentioned in the Introduction, in the past the matrix elements of this oper- 
ator have been computed in the static approximation (the nucleon is assumed 
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to be at rest before and after the interaction with the photon). In the present 
work we go beyond this approximation, since we compute the matrix elements 
using momentum eigenstates. In principle these should be obtained by boost- 



ing |14| the nucleon zero momentum eigenstate ( p0|) . However, the technical 
difficulties associated with boosting are prohibitive. We approximate this op- 
eration by a Peierls-Yoccoz projection, i.e. we consider our model state with 
momentum q to be given by 



\N(q) ) - \f&W¥)^ 



(P q *jM\P q *JM) ' 



(221 



where the square roots are just normalization factors and 

\*J M )=PjM\*Kh). (23) 



The approximation involved in assuming projected instead of boosted states 
is valid for small q. 



Before presenting the formalism to compute the form factors as matrix ele- 
ments of the electromagnetic current ( 21) taken between nucleon states, let 
us recall that 



Jdz [F{z)U{z),R{n)] = Q, 



(24) 



if F(z) is a scalar function of z fL5| . This can be seen writing explicitly the 
commutator, rotating the argument of F and U and redefining z' = TZ(Q)z. 
The previous commutation relation will be very useful in the following. 

Another important point is that, due to the symmetry of the hedgehog, ro- 
tations of this state in spin or isospin space are equivalent. Therefore, as it 
was already pointed out in the previous Section, projecting the hedgehog onto 
spin J implies a simultaneous projection onto isospin T = J, and the two 
projections are equivalent. Hence, the following relations hold: 

PqPjM^hh >= PqPTM^hh >= PTMPq\^>hh >^ Pj M Pq\^hh >, (25) 



where the commutation between operators working in isospin space and oper- 
ators working in ordinary space has been used. The projector Ptm is defined 
similarly to the projector Pjm (Eq. (jl9|)), but with the rotation operator R 
acting in isospin space, replacing the rotation operator in spin space. We shall 
exploit relations ( P5| in the evaluation of the magnetic form factors. 
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4- 1 Electric form factor 



From the definition of the electric form factor (]?]) and using the correspondence 



(22), the electric form-factor is given by 



r , 2 , E fdx(P^ JM \J° m (0)\P-^JM) 

GEiq ] ~ M~ N (P^ JM \P^ JM } ' (26) 



and using the explicit form of the linear momentum projector ([18]) one obtains 
E fdxdbdb'e-^n (U(b')^ JM \J° em (0)\U(b)^ JM ) 



G E (q 2 



M 



N 



f dbdb' e -^- b ')-§ (U(b')*j M \U(b)*j M ) 



(27) 



In principle the form factors should be functions of q 2 only. However, due to 
the approximate treatment of the center-of-mass motion this is in general no 
longer guaranteed and there is a spurious dependence on the angle between 
q and the quantization direction. However, if J = 1/2 states are considered, 
it is possible to show that the form factor is indeed a function of q 2 only. In 
fact, due to parity, the Fourier transform of the matrix element of the current 
has to be a function of (q ■ J) 2 which, for J = 1/2, is proportional to q 2 . We 
can therefore integrate the direction q, both in the current matrix element 
and in the normalization factor at the denominator. After the expansion of 
the exponentials in spherical waves, only the £ = wave contributes and one 
gets 



G E (q 2 ) = 



E J dx dbdb' jo ( 


b+b' 

2 




(U(b')VjM 


J em (0)\U(b)V JM ) 


M N f dbdb' j 


(\^\q) (U(b')*j M \U(b)* JM ) 



.(28) 



The integration on q allows for further simplifications. In fact, we can now 
prove that the translation operator U(b) and the rotation operator R(Q), 
which enters the projector on angular momentum [Eq. (|T9|)1, can be exchanged 
in the previous formula, although they don't commute. Let us, first of all, note 
that 

u\b')J° em (0)U(b) = J° em (b')U(b - b') (29) 



and expand the spherical Bessel functions jo in power series of \b ± b'\q/2. 
Notice that the integration on the direction of the momentum transfer has 
eliminated the dependence of the integrand in Eq. ( P5|) on the angle between 
q and b ±b' . After the introduction of the new variables 

z = b-b' , y = b' (30) 
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we can define the scalar function 



Jdy y 2m+s (cos 9 zy yj° ern (y) (31) 



where 9 zy is the angle between the directions z and y. The l.h.s. of Eq. ( |28| ) 
can be written as a sum of terms of the form: 

J dz (y hh \Pj M F s>Tn (z)U(z)Pj M \y hh ) z 2k+s . (32) 



We can now apply Eq. ( p4|) to each of the previous terms so that one of the 
two projectors on angular momentum in Eq. (p8|) can be eliminated both in 
the numerator and in the denominator. The electric form factor can finally be 
written as 



G E (q 2 ) = i ^ / da J dr J d 3 Q jo^D^fi) 



a ,„ " 



x < U(--)* hh \J» m (r)\U(-)R(n)* hh > (33) 

where the new variables a and r are related to the previous b and b' through 

b' = r-a/2 (34) 

b = r + a/2 (35) 

and the normalization factor reads 

M = J da J d 3 Q Jo (^)V(; M (Q) < * hh \U(a)R(n)\V hh > . (36) 

The numerator of Eq . (|33|) is the sum of three pieces [see Eq. (0)]: the isoscalar 
quark, the isovector quark and the isovector pion contributions. 

4-2 Magnetic form factor 



From the definition of the magnetic form factor (§) and taking again the 
correspondence (|22"D we can write 



. G M {q 2 ) ( x E Jdx(P^j M \J em {Q)\P_^ JM } 

1 2M N [aXq) M N (P^jm\P^jm) [ ' 



where we have introduced a. = (a). 
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In the case of the magnetic form factors, where the space components of the 
electro-magnetic current appear, it is not possible to define a scalar function 
as we did in Eq. fl3~T|). The projection operators will therefore always rotate 
in a non-trivial way the current. To simplify the expression of the magnetic 
form factor we will instead make use of the relations (|25|). The current matrix 
element can therefore be rewritten as 



(* hh \P JM PzJem {0)P-lPjM\*hh) = (*hh\PTMPzJem(0)P-zPTM\y h h) 

= (^ hh \P,_P TM J em (0)P TM P_^ hh ). (38) 

The current contains an isoscalar and an isovector piece. The first one com- 
mutes with the current, while the isovector one transforms as [L6||: 

+1 

PmM JemfiPMM = E ^ q(T M)J em QP m _q M (39) 



where J™ m q stands for the spherical isospin component Q of the isovector 
part of the vector electromagnetic current and 

C Q (T,M) = (10; TM\TM)(1Q; T M—Q\TM). (40) 



To extract the magnetic form factor out of Eq. (|37|) we multiply both terms 
of the equation by a x q and integrate over q. As for the electric form factors, 
this integration is trivial because the Fourier transform of the matrix element 
depends only on q 2 . 

After a straightforward algebraic derivation on obtains the isoscalar (is) part 
of the magnetic form factor: 



x E «w* r 3 [J i : m {r)}k \U(~)R(n)* hh ), (41) 

jk 

where [JgrJfc stands for the cartesian k component of the vector electromag- 
netic current (isoscalar part) and the normalization factor M is given by Eq. 
as before. 



For the isovector (iv) part of the magnetic form factor we get 
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2M N M m n 4r Q 



da / dr 



2gr M ~Q^ ' 



x £«y* ^[Jt, Q (r)] fc |c/(^)^)^} . 



(42) 



5 Results and discussion 



The nucleon electric and magnetic form factors are presented in Figs. 1-4. 
The experimental data shown were taken from Refs. |T7Hl9|. We present our 
results in the space like region < q 2 < 0.5 GeV 2 . 

All the results in Figures 1-4 were obtained using G = 0.2 GeV (actually 
g = 0.024 GeV and M = 1.7 GeV but, as mentioned before, the results 
basically depend on the combination G = y/gM). No parameter was fitted to 
reproduce the experimental form factors. 

For the sake of comparison, in Figs. 1-4 we also present the results for the 
form factors computed in the static approximation, i.e. without performing the 
linear momentum projection. This is the traditional approximation considered 
in previous calculations of the form factors in the framework of soliton models 
(see Ref. || for the linear sigma model and Ref. |J for the double hump version 
of the chiral CDM). 

As it can be seen from Figs. 1-2, the electric form factors are rather satisfactory. 
One has to take into account large incertitudes in the experimental analysis of 
the electric form factor of the neutron, which is obtained from scattering on 
the deuteron and depends hence on the wave function of the latter. The effect 
of the projection on linear momentum is particularly relevant in the proton 
electric form factor. 

The magnetic form factors are less satisfactory. This is probably due to the 
weakness of the spin-spin interaction obtainable in this model, at least working 
within the projected mean-field approximation. 

As it appears from the figures, all the computed form factors underestimate 
the data for large q 2 . It can be interesting to note that, also studying structure 
functions, one sees that in the region of large x = Q 2 /2Mu, where the momen- 
tum carried by the quarks is large, the computed quantities underestimate the 



data |20|]. These problems are probably due to the approximate treatment of 
translational invariance and are therefore not so much related to the specific 
model used in this paper. 

The chiral CDM has now been used to compute many different quantities. 
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We can try to summarize the results. The chiral CDM gives good results if 
problems not involving the spin are considered. This can be seen also from 



the computation of the unpolarized structure functions ^(J. Also the study 
of the transition from nuclear to quark matter within this model seems very 
promising, suggesting a smooth transition between the two phases PI and 



giving interesting results for neutron stars [21 



On the other hand, since the model does not contain enough tensor force, it 
provides poor results for observables which involve the spin. The magnetic 
form factors are therefore not totally satisfactory and the polarized structure 
functions overestimate the data |20| , indicating that most of the spin is carried 



by the quarks because the pion is very weak. It is not yet clear whether 
the weakness of the pionic field is intrinsic to the model and therefore other 
degrees of freedom have to be considered, or stems from the approximations 
used to solve the field equations. Concerning this second possibility, there are 
indications that a large A-N mass splitting could be obtained using the same 
ingredients considered here but allowing the scalar and the vector diquarks to 
have different radii [E3 . 
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Fig. 1 Electric form factor of the neutron as a function of momentum transfer 
in the space like region. The dashed line corresponds to the static approxima- 
tion. The full curve corresponds to the full calculation, i.e. combined linear 
and angular momentum projections. 
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. 2 Electric form factor of the proton. Dashed and full curves as in Fig 
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3 Magnetic form factor of the neutron. Dashed and full curves as in Fig 
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. 4 Magnetic form factor of the proton. Dashed and full curves as in Fig 
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